Fluctuation, dissipation, and thermalization in non-equilibrium 

AdSs black hole geometries 



Simon Caron-Huot 

School of Natural Sciences, Institute for Advanced Study, Princeton, NJ 0854-0, USA^ 

Paul M. Chesler 
Department of Physics, MIT, Cambridge, MA 02139, USA [| 

Derek Teaney 

^ Department of Physics & Astronomy, 

psj Stony Brook University, Stony Brook, NY 11794, USA |^ 

^ (Dated: January 2011) 

PLh Abstract 

•/"J We give a simple recipe for computing dissipation and fluctuations (commutator and anti- 

commutator correlation functions) for non-equilibrium black hole geometries. The recipe formulates 
Hawking radiation as an initial value problem, and is suitable for numerical work. We show how to 
package the fluctuation and dissipation near the event horizon into correlators on the stretched hori- 
O zon. These horizon correlators determine the bulk and boundary field theory correlation functions. 

In addition, the horizon correlators are the components of a horizon effective action which provides 
^ a quantum generalization of the membrane paradigm. In equilibrium, the analysis reproduces pre- 

J> vious results on the Brownian motion of a heavy quark. Out of equilibrium, Wigner transforms of 

commutator and anti-commutator correlation functions obey a fluctuation-dissipation relation at 
O high frequency. 

O 

T — I 

X 



* schuot@ias.edu 

^ pchesler@init.edu 

^ derek.teaney@stonybrook.edu 



1 



I. MOTIVATION 



One of the motivations for studying high temperature gauge theories at strong couphng is 
the striking results from the Relativistic Heavy Ion Colhder and the Large Hadron Colhder 
[H E] . Results on collective flow and the energy loss of energetic probes (in particular heavy 
quarks [3]) indicate that the nuclear size is sufficiently large that macroscopic quantities such 
as temperature, pressure and flow velocity, are useful concepts when characterizing heavy 
ion events. A back of the envelope calculation |3] shows that this would not be possible 
unless the typical relaxation time is of order a thermal wavelength tr ~ ti/T, placing the 
QCD plasma in a strong coupling regime. The AdS/CFT correspondence has led to many 
important insights into the nature of strongly coupled plasmas and energy loss [5j. Of 
particular relevance to this work is the computation of the heavy quark drag and diffusion 
coefficient in A/" = 4 Super Yang Mills (SYM) at large N^. and strong coupling IHHH]- Indeed 
heavy quarks in heavy ion collisions exhibit a strong energy loss and a larger than expected 
elliptic flow, which is qualitatively consistent with a small diffusion coefficient [3]. The actual 
interpretation of the RHIC results is more complicated, since radiative energy loss plays a 
significant (perhaps dominant) role at high momentum where the measurements exist [H]. 
The current results from the Relativistic Heavy Ion Collider on heavy quarks are confusing 
and not generally understood. 

The diffusion of heavy quarks in AdS / CFT is also interesting from the perspective of black 
hole physics. Indeed, the primary goal of this paper is to better understand the thermal 
properties of black holes using the diffusion of heavy quarks in A^ = 4 SYM as a constrained 
theoretical laboratory. 

On the field theory side of the correspondence, the diffusion of heavy quarks is the result 
of a competition between the drag and the noise, which must be precisely balanced if the 
quark is to reach equilibrium. In particular, over a time scale which is long compared to 
medium correlations, but short compared to the equilibration time of a heavy quark, the 
heavy quark is expected to obey a Langevin equation 

I = + (1.1) 

where the drag coefficient r] and the random noise ^ are balanced by the fluctuation- 
dissipation relation 

{C{tW{t')) = 2Tr^5^^{t-t'). (1.2) 

Previously it was shown how this Brownian equation of motion is reproduced by the cor- 
respondence pm [TT] . In AdS/CFT a heavy quark is dual to a long straight string which 
stretches from the boundary to the horizon. At a classical level the straight string is a 
solution to the EOM and doesn't move. This is not the dual of a heavy quark in plasma. 
Fig. [T]^a) shows the geometry of black hole AdS together with a long straight string. In 
our AdS conventions the horizon is at r = 1 while the boundary of AdS is at r = oo. The 
stretched horizon (see below) is at r^^ = 1 + e. Hawking radiation from the horizon causes 
the string to flip-flop back and forth stochastically as exhibited in Fig. [l|^b). The random 
tugs from this flip-flopping string give rise to a random force on the boundary quark which 



is related to the dissipation by the Einstein relation, Eq. (1.2). 

The derivation of this result left much to be desired. In Ref. [lOj it was simply assumed 
that the modes are in equilibrium at the Hawking temperature. With this assumption it is 
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FIG. 1. Two figures from Ref. whicli motivate tliis work, (a) A scliematic of a classical 
string in AdSs corresponding to a heavy quark. The horizon is at r = 1 in the coordinates of this 
work. The stretched horizon is at r/^ = 1 + e and the endpoint of the string is at the boundary 
r-m with Tm JS> 1. Gravity pulls downward in this figure, (b) Hawking radiation from the horizon 
induces stochastic motion of the string in the bulk which we show for three subsequent time steps, 
^1)^2 1^3- The random string configurations give rise to a random force in the boundary theory. 
The Hawking radiation is encoded in an effective action on the stretched horizon r/j = 1 + e. The 
string fluctuations are small, x^y^s ~ 



not difficult to show that the commutator and anti-commutator of string correlations are 
related via a bulk version of the Fluctuation Dissipation Theorem (FDT). When this bulk 
FDT is translated to the boundary theory, the bulk FDT leads to Eq. (1.2). While this 
derivation is physically reasonable, the calculation provides little guidance to out of equilib- 
rium geometries. In Ref. [H], the bulk FDT was derived following a rather complicated and 
un-intuitive formalism [121 [13] • This formalism involves analytically continuing modes as is 
typical in many derivations of Hawking radiation |14f[T6] . The primary goal of this paper 
is to provide a much more physical derivation of the bulk FDT based upon propagating 
near-horizon quantum fluctuations from the distant past up to the bulk. 

In Ref. [H] the effects of Hawking radiation were packaged into a horizon effective action. 
This effective action dictates the dynamics of the fields at r = 1 + e, and provides a quantum 
generalization of membrane paradigm. Although the derivation of the membrane effective 
action involved a complicated analytic continuation, the final form of the effective action 
is very natural. The classical part of the effective action can be determined using the 
classical membrane paradigm, while the quantum part of the effective action is dictated by 
the classical dissipation, and the fluctuation dissipation relation. Once the horizon effective 
action is written down, a short exercise shows how the horizon fluctuation dissipation relation 
leads to the FDT in the bulk and boundary theories. Since the FDT is a direct consequence 
of the fact that the density matrix is exp(— the extent to which this relation holds 
provides an unequivocal measure of equilibrium in the bulk geometry. 

In this paper we will determine the horizon effective action by solving equations of motion 
with appropriate initial conditions rather than analytically continuing modes. Since Fourier 
transforms and analytic continuations are never introduced, it is possible to apply these same 
techniques to non-equilibrium geometries. We also study the thermalization of a string in a 
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non-equilibrium geometry. Furthermore, while we focus on string fluctuations in this paper, 
our analysis easily generalizes to other flelds such as the graviton. 



II. PRELIMINARIES 



In this subsection we will give a brief summary of some of the results of Ref. [TT] in order 
to establish notation. The metric of the black hole AdS space is 



(2.1) 



where the horizon is at r = 1 and the boundary is at r = cxd. L is the AdS radius, 
/(r) = 1 — and T is the Hawking temperature, r is a dimensionless coordinate which 
measures energy in units of temperature. We will also use Eddington-Finkelstein (EF) 
coordinates to describe the near horizon dynamics. In this coordinate system the metric is 
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A{r)dv'^ H — —drdv + r^dx^ 
ttT 



where A = r"^ f{r) and v is the EF time 



t + 



ttT 



dr 

J^2 



(2.2) 



(2.3) 



Ingoing lightlike radial geodesies have v = const, while outgoing lightlike radial geodesies 
satisfy dr/dv = 'irTA/2. From now on we will set the AdS radius to one, L = 1. 

For simplicity, consider fluctuations along an inflnitely long straight string {i.e. an in- 
finitely massive quark) stretching from the horizon to the boundary. The stationary bound- 
ary endpoint is at x = 0, and small fluctuations around the straight string solution are 
parameterized by x{t,r), where x denotes displacement of the string in the x direction. Ei- 
ther t and r or v and r are taken to be the world sheet parameters, t = and r = a^. The 
action of these world sheet fluctuations is derived by linearizing the Nambu-Goto action: 



S 



27 



dtdr 



-IVhh'^'duxd, 



X 



(2.4) 



where /i, u run over t, r ot v,r depending on the coordinate system. For example, the world 
sheet metric in EF coordinates is 



—A{r)dv^ H — — drdf 
ttT 



(2.5) 



We note that the drag coefficient of the heavy quark (see Eq. (1.1)) is related to the coupling 
between the metric and world sheet fluctuations [6HH] 
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The analyses in the following sections make this physical interpretation of 7] clear. 



(2.6) 
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The goal of this paper is to show by solving equations of motion that in equilibrium the 
retarded propagator, 

iGraitin\t2r2) = 0{t - t') {[x{ti,n),x{t2,r2)]) , (2.7) 
and the symmetrized propagator, 

Grr{tiri\t2r2) = ^ {{x{ti, n), X{t2, Ts)}) , (2.8) 

are related by the fluctuation dissipation theorem 

Grr{uJ,r,f) = -(1 + 2n(w))ImG',„(w,r,f) . (2.9) 

Here n{uj) = l/{exp{uj/T) — 1) is the Bosc- Einstein distribution function. This relation is a 
direct consequence of the fact that the density matrix is exp{—H/T) and signifies that the 
fluctuations are in equilibrium with the black hole at temperature T. 

The advanced propagator is related to the retarded propagator by time reversal, 

iGar{tiri\t2r2) = iGra{t2r2\tiri) , (2.10) 

while the spectral density is the full commutator 

iGra-ar(tiri\t2r2) = ( [£(tl , Ti) , x(t2, ^2)] ) . (2.11) 

Simple manipulations show that 

- 2lmGra{i^, r, f) = iGra-av{^-, T, f) , (2.12) 

and thus the fluctuation dissipation theorem is a relation between the commutator and anti- 
commutator which signifles equilibrium. Both the commutator and the anti-commutator 
will be determined by solving equations of motion with appropriate initial conditions. 



III. EQUILIBRIUM STRING FLUCTUATIONS IN AdS^ FROM EQUATIONS OF 
MOTION 

A. Equations of motion and boundary conditions 

Let us analyze the equations of motion in more detail. Recall that the retarded correlator 
is a Green function of the equations of motion 



y5 

27r 



df.g.,Vhh^''dJ\ Gra{tiri\t2r2) = S{h - t2)S{n - rs) , (3.1) 



and is required to vanish when tiVi is outside the future light cone of t2'"2- By contrast, 
the full commutator {i.e. the spectral density) is not a Green function but satisfies the 
homogeneous equations of motion 
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d;, V^/l'^"^,] Gra-raitm \t2r2) = , (3.2) 
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where the initial conditions are determined by the canonical commutation relations. Simi- 
larly, the symmetrized correlation function also satisfies the homogeneous equations of mo- 
tion 



2tx 



df^g^xVhh'^'d^ Grritiri\t2r2) = 0, 



(3.3) 



but the initial conditions are determined by the density matrix of the quantum system 
far in the past. The appropriate initial conditions for Grr and Gra-ar are discussed more 



fully in Section III C Finally all bulk to bulk correlation functions {Gra, Gr 



satisfy 



Dirichlet, or normalizable, boundary conditions for asymptotically large radius, i.e. G — )■ 
for ri, r2 — oo. 

Since the supergravity equations of motion are essentially coupled oscillators, it is useful 
to recognize that the retarded propagator for the simple harmonic oscillator is independent 
of the density matrix. Only symmetrized correlations depend on the density matrix and 
reveal a thermal state. Since the simple harmonic oscillator clearly illustrates the role of 
the density matrix, we show how to compute commutator and anti-commutator oscillator 
correlations using the Keldysh formalism in Appendix |Xj 



B. Horizon Correlators 



The equations of motion propagate initial data in the past to the future. This can be 



made manifest for the symmetrized correlator by writing down a formal solution to Eq. (3.3) 
in terms of retarded correlators. Specifically, given Grr and its time derivatives on some time 
slice ti = t2 = to, the solution to Eq. (3.3) at later times is given by 



ar(i|2) 



2n 



J dr[g^^Vhh'\r[)Grail\l')^, 



X 



2n 



dr', g,^Vhh'\r',)Grai2\2')W 



Grr{l'\2'), (3.4) 



where t[ and t2 are set equal to to after differentiating, and d = d 



This formula 

expresses the uniqueness of the correlator given its value and time derivatives on a Cauchy 
surfacej^ The physical interpretation of this solution is easy to understand. The two retarded 
Green functions appearing in the integrals are convoluted with the separate arguments of 
the initial data. These retarded Green functions causally propagate the initial data forward 
in time. 

To gain qualitative insight into how initial data is propagated by the retarded Green 
functions. Figure [2] shows a congruence of outgoing radial null geodesies starting at time t = 
to- A generic geodesic reaches the boundary in a time At ~ 1/T. The information which is 



The formula manifestly satisfies the equations of motion (3.3) for ti,t2 > <o- To see that it satisfies the 



boundary conditions in the limit ti — > tQ, one must know the time derivatives of Gra(l|l') for ti — > t[. 
This derivative can be obtained by using the fact that G^a (1|1') vanishes for ti < t'l, and by integrating 
ti across t'l with the equations of motion ( |3.1[ ) to yield the canonical commutation relations 

hm ^g^^Vhh'^dtMM'^') = S{n - r[) . 
tl—^t[ Zn 

An analogous formula holds for the derivative with respect to t'l. 
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Relevant initial data 
at late times 
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r = 1 + e 
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FIG. 2. A congruence of outgoing null radial geodesies starting at time to- Qualitative insight on 
the propagation of initial data specified on the slice at t = to can be understood from the congru- 
ence. Generic geodesies reach the boundary in a time At ~ 1/T. Initial data propagated on such 
trajectories reflects from the boundary and falls into the horizon with an infall time of order 1/T. 
Geodesies originating exponentially close to the horizon take much longer to escape. Consequently, 
at late times the above-horizon geometry is filled with geodesies emanating exponentially close to 
the horizon at t = to- Because of this, the only initial data relevant at late times consists of the 
initial data exponentially close to the horizon at t = to. 



propagated along such trajectories reflects off the boundary and falls into the black hole with 
an infall time also of order 1/T. Thus, at times At ^ the only outgoing geodesies which 
populate the geometry above the stretched horizon are those which started exponentially 
close to the horizon at t = to. Moreover, the initial data propagated from this exponentially 
narrow strip to the above-horizon geometry will be dramatically redshifted. Because of the 
redshift, the only finite wavelength contributions to the symmetrized correlator will come 
from the UV part of the initial data near the horizon. This UV part simply comes from 
coincident point singularities in Grr, which encode quantum fiuctuations in the exponentially 
narrow strip near the event horizon, r = 1. 

This discussion suggests a strategy for computing the symmetrized correlator at times 
long after to. First, all initial data above the stretched horizon at r^j = 1 + e can be neglected 
as it gets absorbed by the black hole at times t ~ to + 1/T. Initial data which lies in the 
strip between the event horizon and the stretched horizon can be evolved until it reaches 
the stretched horizon, where it determines an effective horizon correlator. This correlator 
encodes quantum ffuctuations in the past and sources radiation on the stretched horizon 
that subsequently propagates up to the boundary. This is depicted graphically below in 
Figure |3} Given the discussion of the preceding paragraph, we expect that the form of the 
horizon correlator will be independent of details of the initial data specified in the distant 
past. 

The effective horizon correlation function can be found by exploiting the composition law 
obeyed by retarded propagators. Let S'ra (^i'"2 1 ^2'"2) denote the retarded propagator in the 
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region 1 < r < 1 + e subject to the reflective Dirichlet conditioE^ at r = 1 + e 

graihri = rh\t2r2) = 0, 5'm(^l'^l|^2^2 =rh) = 0. 



(3.5) 



Appendix |B] analyzes the Wronskian of the retarded propagators Gra and Qra, and determines 
the following composition rule 



GraiMl') 



dt2 Grail\2) 



V 



Vhh''-{r2)dr, 



9ra{2\l'), 



(3.6) 



r2=rh 



where tiVi is outside the stretched horizon, while t[r[ is inside the stretched horizon. This 
identity, schematically depicted in Fig. [3} is the mathematical statement of how information 
is propagated up from near the horizon to the stretched horizon, and then up from the 
stretched horizon towards the boundary. Substituting this composition law into the solution 
(3.4) and neglecting contributions to the integrals from r > 1 + e, we find that above the 
stretched horizon Grr takes the form 



Grr{tiri\t2r2) = / dt;dt'2 [~Gra{tiri\t[rh)] [-Gra{t2r2\t'2rh)] G%{t[\t'2 



(3.7) 



where the "horizon correlator" G^^ is determined only by the dynamics between the horizon 
and the stretched horizon: 



Grriti\t2 



r\,T2=rh 



(3.8) 



Here grr{ti^2\t2f2) is the solution to the homogeneous equations of motion (3.3) with reflec- 
tive Dirichlet boundary conditions at the stretched horizon, together with the prescribed 
initial conditions close to the horizon at t = tojj We will always denote bulk correlation 
functions inside the strip with lower case letters^] 

While the symmetrized correlator measures the degree of occupation of microstates, the 
spectral density iGra-ar{<^) measures the density of available states. Since the spectral 
density is also a solution to the homogeneous equations of motion, the discussion of the 
previous paragraph can be repeated mutatis, mutandis yielding 



iGr 



(tiriltars) = / K^t'^ [-Graitin\t\rh)] [-Gra{t2r2\t'2rh)] ^G'^.^_^^.{t\\t 



21 1 



(3.9) 



where the horizon spectral density G^a-ar is related to the spectral density inside strip 
as in Eq. (3.8) 



-r]y/hh"\r{)dA \-r]y/hh"\r2)dr 



(3.10) 



gra-aritiri\t2r2) 

In this equation gra-ar is subject to the same reflective Dirichlet boundary conditions at 
r = 1 + e. We will see that the horizon correlators, G^^_^^ and G^^, form the components of 



a horizon effective action, which we will study more completely in Section IV 



^ We note that it is not necessary that gra satisfy reflective Dirichlet condition at r = rh- This choice is 

simply a matter of latter convenience. 
^ In particular, (7rr(l|2) is given by the same expression in Eq. (3.4|, but with the replacements G',,.(l|2) — )• 

5(rr(l|2), Gra(l|2) — gra(l|2) and with the limits of integration running from 1 < r < 1 + e. 
^ This correlator will always be written with its arguments grr{viri\v2r2), and can not be confused with 
the metric coefficient grrifjv). 
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to t > 



FIG. 3. The composition law for retarded Green functions. The stretched horizon separates the 
two steps of the evolution. From the perspective of an observer in the exterior, the strip between 
the horizon at r = 1 and the stretched horizon at r = 1 + e produces a simple horizon effective 
correlator, which acts like a source of radiation from the stretched horizon. 



Eq. (3.7) and Eq. (3.9) are the key equations in our study of thermalization. Together 
they show that emission and absorption of fluctuations in the bulk can be encapsulated into 
a horizon fluctuation and a horizon resistance. These equations show that we can focus our 
attention on the dynamics inside a small strip between the horizon and stretched horizon. 
Thermalization inside the strip will be simple to understand analytically. 



C. Propagating correlations to the stretched horizon 



1. The horizon symmetrized correlation 



As is evident from Eqs. ( |3.7 ) and (3.8), to determine the symmetrized correlator we must 
flrst determine grr{tiri\t2r2) in the strip between the horizon and the stretched horizon. For 
describing the evolution inside the strip it is useful to use Eddington-Finkelstein coordinates. 
In these coordinates the equations of motion for r < rh are 



and 



dri 

d_ 
dr2 



2h^ 



dvi 



d 



OV2 



_d_ 
dri 

_d_ 
dr2 



grr{viri\v2r2) 



grr{viri\v2r2) 



0, 



(3.11) 



(3.12) 



where the metric coefficients are given in Eq. (2.5). Since we are interested in the ultra-violet 



irregular solution to Eq. (3.11) and Eq. (3.12), we have neglected the radial derivatives of 
(yfa-a,, which are small compared to the radial derivatives of (^rr ('yi'"ib2^2)j3 With this short- 

^ A posteriori one can verify that the neglected derivatives, e.g. (drgxx) {dvgrriviri\v2r2)), are small com- 
pared to the terms which are kept, e.g. gxx drK"^ dygrr^'^) ■ 
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distance approximation, the linear operator in Eq. (3.11 ) becomes a product of two operators, 
dr and + ^Adr- Recall that in Eddington-Finkelstein coordinates, outgoing null radial 
geodesies satisfy ^ = ^A, whereas infalling null geodesies satisfy v = const. Thus, the 
derivatives act along these geodesies, implying that functions of e~^'^'^^(r — 1) and v satisfy 
the equations of motion in the near horizon limit. The general solution is therefore 

grrivmM = fiie-'-'^^^in - i),e-2-^^Hr2 - 1)) + /2(e-2-^^Hn - l),V2) (3.13) 

+ /3(^i,e-^'^^^Hr2-l)) + /4(t;i,t;2), 
where the fn{xi,X2) are arbitrary functions. Requiring that the boundary conditions 

9rr{vi,ri = rh\v2r2) = 0, and grr{viri\v2,r2 = rh) = , 
are satisfied at all times, we find that 

grr{viri\v2r2) = + 

-fie 
-fie 

+ fie'''^'''''irH - l),e-^-^^'^Hr, - 1)) , (3.14) 

where f{xi,X2) is determined by initial conditions. 

The above solution is a linear combination of modes which are outgoing from the hori- 
zon and modes which are infalling towards the horizon. The infalling modes are in fact a 
consequence of the outgoing modes, as the reflective Dirichlet boundary conditions at the 
stretched horizon turn outgoing modes into infalling modes, which are subsequently ab- 
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(ri - 
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-2itTv2 
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~2ttTvi 


(ri - 
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e~ 


-2ttTv2 


[rk. 


-1)) 


~2ttTvi 


i^h - 


1) 
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-2nTv2 


ir2 


-1)) 




irh - 


1) 
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'2ttTv2 


irh 


-1)), 



sorbed by the black hole. Inspection of Eq. (3.14) shows that if f{xi,X2) is analytic, this 



reflection and absorption would always lead to grri'Viri\v2r2) — )■ at late times. Further- 
more, all outgoing modes originate exponentially close to the horizon. This follows from the 
fact that as e~^'^^^(r — 1) is constant on outgoing null geodesies, in the limit v — >■ —oo we 
must have r — )■ 1. Therefore, for evolution near v = 0, relevant initial data specified in the 
distant past will come from an exponentially narrow strip which is exponentially close to 
the horizon at r = 1. In this exponentially narrow strip grriviri\v2r2) is not analytic and 
contains a logarithmic coincident point singularity 

grriviri\v2r2) = log |/i(fi — f2)(?"i — ^2)1 + exponentially small terms, (3.15) 

where is a constant]^ 

By matching Eq. ( |3.14[ ) on to Eq. (3.15) we conclude that we must have 



/(xi, X2) = log |xi — X2I + exponentially small terms. (3.16) 



^ The coincident point limit can be found be computing the anti-commutator for a free massless scalar field 
theory in l+ld flat space 

^ mx),m}) = log if^v^.x^x-'i . 



Eq. (3.151 is found by the principle of equivalence and by comparing the normalization of the string action 



(2.4 1 to the canonical normalization. 
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Substituting Eq. (3.16) into Eq. (3.14), we find that long after initial conditions were specified 
and up to exponentially small corrections, we have 



9rr[Viri\V2r2) = - 



+ 



log \e-'-^^^{r,-l) - (r2-l)| + log je" 



■2-itTAvi 



ri-1) - (rh-1) 



log |e 



-2TrTAv 



log|e-^'^^^^(r,-l)-(r,-l)|. 

(3.17) 



In writing Eq. (3.17) we have pulled out a common factor of e '^'^'^^'^ which cancels between 



the different terms due to the reflective boundary conditions. 

The fact that grr{viri\v2r2) is time translationally invariant at late times is due to the 
competition between emission and absorption. Quantum fluctuations give rise to the loga- 
rithmic coincident point singularity in grr{viri\v2r2), and this singularity acts as a constant 
source of radiation. Thus, radiation is continually produced, redshifted, reflected off the 
stretched horizon, and absorbed by the black hole. This dynamical process reaches a steady 
state which determines grr{viri\v2r2) at late times. As we will discuss in Section HID, this 
steady state is in fact a thermal state. 

With grr{viri\v2r2) known, we may compute the horizon correlator G^^{vi,V2) via 



Eq. (3.8). The result reads 



G'^,{v,,V2)=-'^d,A2 log|l-e-2-^(^^-^^)| 



(3.18) 



2. The horizon spectral density 

As in the case of grr{'ti^i\t2f2) , it is convenient to use EF coordinates [v, r) to determine 
gra-ar(tiri\t2r2) ■ The spcctral density gra-ar{viri\v2r2) obeys the same equation of motion 
as grr{viri\v2r2) and has the same reflective boundary conditions at the stretched horizon. 
Thus, the general solution for gra-ar{viri\v2r2) is similar to the corresponding solution for 



S'rr (^1^1 1'y2'"2) Eq- (3.14). Explicitly, we have 
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{r-h - 


1), 


e' 




-1)), 



gra-ar[Viri\V2,r2} 



(3.19) 

where fra-ar is determined by initial conditions. 

In contrast to the symmetrized correlator, where the initial conditions are determined 
by the density matrix in the past, the initial conditions for the spectral density are state 
independent and are determined by the canonical commutation relations. Thus, it is not 
necessary to evolve for a long time before reaching a steady state solution. 
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For viTi close to f flat-space physics determines gra-ar 

gra-ar{viri\v2r2) " 7" (sign(fi - V2) - sign(ri - r2)) 

Arj 



(3.20) 



Comparing Eq. (3.20) to Eq. (3.19), we conclude that for any vi and V2 we must have 

1 



ra—ar V^i ) 



Xi,X2 



4?7 



signi^xi - X2j 



(3.21) 



The infalling and outgoing modes cancel at spacelike separation due to the sign function. 
One can verify (with careful algebra) that the Qra-ar satisfies the canonical commutation 
relation 

r]Vhh^\ri) lim dt^gra-ar{tiri\t2r2) = 5{ri - ra) . (3.22) 

t2— >ti 



Substituting gra-ar into Eq. ( |3.10[ ), we determine the horizon spectral density 

Gt-ar{Vl,V2)=2r]6'{Vi-V2). (3.23) 



D. The Bulk Fluctuation Dissipation Theorem 

The interpretation of the above results becomes clear in Fourier space. The Fourier 
transform of the horizon symmetrized correlator is 



Gt(co) 



dv e*^^ a" log 11- e- 



-2ttTv I 



r]u (1 + 2n(w)) 



(3.24) 



where n{u) = l/(exp(a;/T) — 1) is the Bose-Einstein distribution. The Fourier transform of 
the horizon spectral density is 

tGl_,,{u) = 2r]u. (3.25) 

Thus, after the decays of transients in the initial data, the horizon correlation functions obey 
the fluctuation dissipation relation 

G%{u) = iGl_,,{u) Q + n{uj)^ . (3.26) 

Now we can see how the bulk will thermalize from the se h orizon correlations [Uj. In 
Fourier space the expression for Grr in terms of Gj?^ in Eq. (3.7) becomes 



Grr{u,ri,r2) =Gra{(^,ri,rh)Grai-u},r2,rh) [G%{uj)] , 
and the spectral density obeys a similar equation 

iGra-ar{^^,ri,r2) = Gra{uj , n, rh)Gra{-^^ , r2, Th) [iG'^a-ari^)] ■ 



(3.27) 



(3.28) 



^ The coincident point limit can be found by computing the commutator for a free massless scalar field 
theory in 1+ld flat space 

mX), 0(0)]) = -h{-Tj^,Xf^Xn sign(t) = ^(sign(t + z) - sign(-t + z)) . 

The sign function is the only functional form which allows left-moving and right-moving modes to can- 



cel at spacelike separations. Eq. (3.20) is found by the principle of equivalence and by comparing the 



normalization of the string action (2.4) to the canonical normalization. 
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FIG. 4. Feynman graph used for computing the symmetrized correlation function ri, r2), 



see Eq. (3.27). 



Thus the horizon fluctuation dissipation relation trivially implies the relation in bulk 



G'„(a;,ri,r2) = iGra-ari(^,ri,r2) [ - + n{u) 



(3.29) 



In Section IV we will show how the convolution in Eq. (3.27) is the result of coupling the 



horizon effective action to the bulk. Anticipating these results, Fig. ffl shows the Feynman 



graph corresponding to Eq. (3.27). 



Physically, the bulk fluctuation dissipation theorem follows from its horizon counterpart 
because any fluctuation in the bulk must have crossed the stretched horizon at some point 
in the past. Previously the form of the horizon correlators was derived either by assuming 
equilibrium, or by using a complex set of analytic continuations. We see that it is the 
endpoint of a simple competitive dynamics inside the strip. 



IV. THE HORIZON EFFECTIVE ACTION 



It is enlightening and useful to package the above manipulations into a path integral 
formalism. The claim is that these steps correspond to integrating out the fields inside the 
stretched horizon. The horizon correlation functions defined above are the components of 
the resulting effective action. 



A. The Keldysh Formalism 

In any quantum statistical system, correlations are computed by tracing the density 
matrix p over the Heisenberg operators. This is the entire purpose of the Keldysh contour 
[T7j . The generating function of string correlation functions in bulk is 



Z[J-i,J-2] = Tr 



X2{to,r) 



[Dxi][Da;2; 



gi5i-i52 gi / dMrJ-i(t,r)xi(t,r) ^- 



■ifdt'dr'T2{t',r')x2{t'y) 



(4.1) 
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xi(t,r) 

to ) ^max 

X2(t,r) 

FIG. 5. The Keldysh contour from an initial time to where the density matrix is specified to a 
final time, tmax- Xa{t, r) = xi{t, r) — X2{t, r) must be zero at tmax- 

where [Dxi] indicates a bulk path integral, i. e. 

[Bxi] =Y[dxi{t,r) . (4.2) 

t,r 

The path integral is defined along the Schwinger-Keldysh contour shown in Fig. [5} where the 
"1" type path integral is the amplitude of the process, while the "2" type path integral is 
the conjugate amplitude of the process. The trace is over the initial density matrix p which 
determines the initial conditions xi(to,r) and X2ito,r). T\ and T2 are sources, which in this 
case are simply the external forces applied to the string. Variation of generating function 
with respect to T\ and T2 yields time ordered, anti-time ordered, and Wightman correlation 
functions [Hj. 

Instead of using the 12 variables, we will rely on a rewritten version of the Keldysh 
formalism known as the ra formalism which is dramatically simpler. We define retarded (r) 
and advanced (a) fields and sources 

_X\\ X2 _ Tr_*^l+*^2 _ , . 

Xy — ^ , Xa — Xi — X2 , J- r — ^ 7 J'a — J'\~J'2i l^-'Jj 

Since Xa encodes the the differences between the amplitude and conjugate amplitudes, Xa is 
a small parameter in the classical limit [18]. In terms of r and a fields the action becomes 



5i - ^2 = ^ / dtdr \-^fhh^''d^Xrd,x^ . (4.4) 

The two point functions in the "ra" formalism are familiar and explain the notation of 
Section lU 

iGra{t\r\\t2r2) = {Xr{ti,ri)xa{t2,r2)) = 0{ti - t2) {[x(ti,ri),x{t2,r2)]) , (4.5) 

Grr{tiri\t2r2) = {Xr{ti,ri)Xr{t2,r2)) = ^ ( {x(ti , Ti) , X (^2 , ^2) }) . (4.6) 

The causal structure of quantum field theory is rendered transparent in ra formalism. 
Since at tmax along the contour, Xi(tmax) = 2^2 (^max), the path integral must be solved 
with the boundary condition, Xa{t) — )■ for t — )■ 00. Thus, whenever an "a" type field 
is evaluated at a later time than all other field insertions the correlator vanishes. For 
example, Gaa(t,t') = since an a field is always evaluated last. Similarly the retarded 
correlator, Gra{t,t'), vanishes whenever the a field is evaluated at a later time than the 
r field. Gra determines (minus) the retarded linear response to a classical force. We will 
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exhibit the retarded correlator with an arrow to indicate the direction of time recorded by 
this propagator 

^ = iGraiuj,r,f) . 

While the retarded correlators reflect the response to a classical force, symmetrized correla- 
tors encode the fluctuations in the system. Since the symmetrized correlation function does 
not represent causal response to a classical force, but rather a time dependent correlation 
which arose from a specified initial condition, we will notate this correlation (as in Ref. [19j ) 
with 

^ H ^ =G„(a;,r,f) . (4.7) 

Higher point correlation functions involving r and a indices have a similarly simple inter- 
pretation [T7] . 



B. The effective action from the path integral 



To obtain the horizon effective action we integrate out all field fluctuations inside the 



stretched horizon at = 1 + e. The path integral in Eq. (4.1) becomes 



(4.^ 



r>rh 



where the horizon effective action is the path integral 



r<r^ 



iSi-iS2 



(4.9) 



with fixed boundary values on the stretched horizon x'^{t). Here the "s" sub-label denotes 
the Schwinger-Keldysh index. The horizon effective action must be traced over the initial 
data extremely close to the event horizon. The fixed boundary conditions couple the field 
fluctuations with r > 1 + e to field fluctuations with r < 1 + e. The trace determines how 
vacuum fluctuations in the past influence the future dynamics. It should be understood 



that Eq. (4.8) (where the initial density matrix is retained only inside the stretched horizon) 
is an approximation that is valid after most of the initial data has fallen into the hole as 
discussed in Section UTTl 

To quadratic order the effective action can be expanded as: 



dtdix':,{t)[zGra{t,i)]x'^{i)-- / dtdix':{t)[G%{t,i)]x':{t), (4.10) 



where we can determine the horizon correlators via differentiation, e.g. 



1 



i 6S, 



cS 



z2 6x'^{t)6x'^,{t) 
1 iSScs 



^2 (5x^(t) 5x^(t) 



(4.11) 
(4.12) 
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A short calculatior]|^ shows that every derivative —ibjbxa or —ibjbxr brings down a factor 
of —rjy/hh'^^drXr or —rjy/hh^^drXa respectively. Thus 



Gl{t,t)= lim \-r]y/hh''^{r)dr] \-vy^h"'ir)dr \ gAhnlt^r^) 

r,r^rh L J L J 

G^,(t,t)= lim \~r]Vhh'-'{r)dr] \-vVhh^'{f)dr \ gra{tiri\t2r2) 



r,r^rh 



(4.13) 
(4.14) 



where the lower case correlation functions are defined from the path integral for r < rh with 
vanishing boundary conditions at r = r/^ 



(•••); 



Tr 



PI Bxs5[xsiv,r„)]e''''-"'' 



(4.15) 



These relations are familiar from the context of the AdS/CFT correspondence where one 
takes radial derivatives of the fields as r — )■ oo |20l422] . The functions grr{tiri\t2r2) and 
gra{tiri\t2r2) (or equivalently gra-ar = gra — gar) are the same as in section III and their 



relation to G^^ and is identical to that derived in Eq. (3.8). Further the reflective 



boundary conditions imposed on grr and gra appear naturally in the path integral formalism. 
We conclude that the horizon correlators Gj?^ and G^^_^^ defined and computed in the 
previous section are precisely the components of a well-defined effective action. Given a 
procedure for regularizing ultraviolet divergences in gravity, this effective action could be 
computed to any desired order in perturbation theory following the methods of Ref. |23j . 



C. Summary of the horizon effective action 

The effective action in equilibrium was obtained previously by Son and Teaney |llj by 
analytically continuing modes across the horizon. Its form is simple, 



h/ \ 



(4.16) 



where the retarded horizon correlation function is simply 

G'^, = -tuV, (4.17) 
and the symmetrized part of the horizon action obeys a fluctuation dissipation relation 

G'^^ = -(1 + 2n{Lo)) IraG'^^ico) . (4.18) 
Fourier transforming back to time, the horizon action reads 

^ Xait') . (4.19) 



iSes = -% I dtx2{t)r]dtxAt) - ^ / dtdt' x^^it) 



TT 



This is because the only dependence on comes through the boundary terms of the kinetic term, e.g. 
- / dM. ,V-hh-d...d... ^ -,V-hh-x^d... + / dM. x.a. UVhh-dr: 

J r=rh J V 
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This action compactly summarizes all the correlations that appear through quadratic order 
in Xn 



and agrees with the results of Section III D 



The horizon effective action is useful. For instance, since the bulk action in Eq. (4.4) has 



no rr and no aa components, the first Feynman graph which contributes to the correlation 
{xr{vi^ri) Xr{v2-,r2)) ^ is shown in Fig. |4]and is written in full in Eq. (3.27), with the lines 
explained in Section IV A Thus, perturbation theory with the bulk and horizon actions 



(3.9) 



transparently produces the convolution formulas given in Eq. (3.7) and Eq. 

The retarded horizon propagator in the action, G^^, reflects the resistance on the horizon, 
— r^x, and is valid at all frequencies [TT]. Indeed, the classical dissipation encoded by the ra 
part of the action can be derived simply from the classical membrane paradigm. Variation 
of the effective action gives the horizon force 



6S, 



cff 



6xn 



-r]dtXr 



(4.20) 



where the force is J-"/^ = —•q^/hh'^'^dj.Xr as required by the membrane paradigm [24j. In the 
zero frequency limit the force J> is independent of radius [HI EH], implying that horizon drag 
coefficient rj is the same as boundary drag coefficient in the Langevin equation [B]. 



V. NON-EQUILIBRIUM CORRELATORS 
A. Setup 

In this section we wish to show how to generalize the non-equilibrium horizon effective 
action. Interesting non-equilibrium geometries to consider can be found in [261 - I28] . However, 
for definiteness we will consider the non-equilibrium geometry discussed in Ref . [2S] • In this 
work, an excited state in the boundary quantum field theory was created by briefly turning 
on a time-dependent Ad gravitational field, which was taken to be translationally invariant. 
The gravitational field did work on the quantum system, producing an excited state which 
subsequently thermalized. In the dual bd gravitational system, turning on a Ad gravitational 
field corresponds to deforming the Ad boundary of the bd geometry. Before the deformation 
was turned on, the bd geometry was taken to be AdS^, which is dual to the vacuum state. 
The deformation in the boundary geometry produced gravitational radiation which fell into 
the bulk. This infalling radiation resulted in the process of gravitational collapse, changing 
the initial AdS^ geometry to one which had a black hole, and the relaxation of the black 
hole to equilibrium encoded the thermalization of the expectation value of the stress tensor 
in the dual quantum theory. 

Translation invariance allows the bd metric to be written 

ds^ = Adv'^ + 2drdv + {e^dxl + e'^^dxl) , (5.1) 

where all coefficients A,J],B are functions of radius r and time v. The metric coefficient 
A{v,r) together with the outgoing radial geodesies calculated in this geometry are shown in 
Fig. |6]and is reproduced from [215]. Outgoing light like geodesies satisfy dr/dv = A/2. We 
will determine the non-equilibrium string correlators in this transient geometry. 

A salient feature of the outgoing geodesies is their ultimate bifurcation at the event 
horizon, ro{v), which is shown by a thick black line in Fig. [61 This bifurcation is reminiscent 
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FIG. 6. Figure from Ref. ^26j. The congruence of outgoing radial null geodesies. The surface 
coloring displays A{v,r)/r'^ . The excised region is beyond the apparent horizon, which is shown 
by the dashed green line. The geodesic shown as a solid black line is the event horizon; it separates 
geodesies which escape to the boundary from those which cannot escape. 



of Fig. [2j It is convenient to switch coordinates to a system of coordinates where this black 
line is flat 

p =r — ro{v). (5.2) 

Note ro(t>) deflnes a light like outgoing radial geodesic, so that the metric is 

ds'^ = -{A- Ao)dv^ + 2dpdv + T? {e^dx\ + e-'^^dx{) , (5.3) 

where Ao{v) = A[ro{v),v). In the near horizon geometry, where r ~ ro{v) we can approxi- 
mate 

^/ N A r \ dA{r,v) 
A{r,v) - Ao{v) ~ ^ ^ 



dr 



(5.4) 



r=rh(v) 



where here and below we will deflne the stretched horizon at ph = e, i.e. rh{v) = ro{v) + e. 
For future convenience we deflne the "effective temperature"!^ 



27rTeff(w) 



1 dA{r,v) 



dr 



(5.5) 



r=rhiv) 



The action of the string fluctuations in the x direction is the same as Eq. (2.4), but the 
metric coefficients depend on time and radius. With the goal of determining the horizon 
effective action it is useful to deflne a non-equilibrium drag coefficient 



(5.6) 



^ We note that T^ff should not be interpreted as a temperature at all times. In particular, due to the 
teleological nature of event horizons, Tes is non-zero even before the boundary geometry has changed - 
i.e. when the dual boundary quantum theory is still in the vacuum state. It is only at late times when 
the black hole starts to thermalize that Teg can be interpreted as a temperature. 
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B. Calculation 



The computational procedure for computing correlators in the non-equilibrium case is 



remarkably similar to the equilibrium case discussed in Sections III and |IV[ For both Grr 
and Gra-ar oue can write down a solution to the equations of motion in terms of retarded 



Green functions convoluted with initial data, as in Eq. (3.4). Furthermore, for initial data 
specified suitably far in the past, the relevant initial data for evolution near f = (when the 
boundary geometry is changing) will come from a narrow strip near the horizon One can 
then repeat the analysis of Section III B and conclude that Grr and Gm-nr are determined 



by horizon correlators G^^ and as in Eqs. (3.7) and (3.9) 

In the non-equilibrium case G^^ is given by 



Gt{v,,V2) 



-ri{v^)VhhPP{p{)dp} \-v{v2)^hPP{p2)dJ grr{vipi\v2P2) (5.7) 



where grr{vipi\v2P2) is a solution of the homogeneous equations, but is confined to the strip, 
< p < e. Furthermore grriyipi\v2P2) should satisfy the boundary conditions grriyi, pi = 

(^\V2,P2) = grr{ViPi\v2P2 = e) = 0. 

The equations of motion for grr{vipi\v2P2) are 



dpi 



grr{Vlpl\v2p2) = ^ , (5.8) 



OVi Opi Opi OVi opi 

where all the metric coefficients depend on v and r. Without approximation, we have 



d_ 
dpi 



OVi OVi opi 



d{g,,Vhhn d ^d{g^^VhhP- 



dpi 



dvi 



grr{VlPl\v2P2) 

grr{VlPl\v2P2) = ^ ■ (5.9) 



dpidv 



The above equation of motion contains both first and second derivatives. However, as 
we are interested in solutions which are irregular near the horizon, all terms with single 
derivative operators acting on gTr{vipi\v2P2) can be neglected. This approximation leads to 



d_ 
dpi 



2r]\/hk 



, d di^qVhhP") 
dvi 



dvi 



d_ 
dpi 



grr{viPl\v2P2) = . 



(5.10) 



and an analogous equation for V2 and p2- Inspecting the above equation, we se e that 
\/ri{vi)ri{v2)grrivipi\v2P2) is annihilated by the operator dp [d^ + ^Adp] . As in Section III C 1 



the consequence of this is that any function which is constant on null radial geodesies satis- 
fies the equations of motion. Near the horizon this translates to ^/ri{vi)ri{v2)grri'Vipi\v2P2) 
being a function of v and 



Poutip,v) = pe~ 



(5.11) 



To make this more precise, suppose instead of starting off with an initial geometry which was AdS^, the 
initial geometry consisted of a static black hole geometry at temperature Tinitiai = S. Assuming initial 
data is specified at times w ^ —1/(5 in the past, all relevant initial data for future evolution around v = 
will come from an exponentially narrow strip which is exponentially close to the horizon. Of course, one 
can always consider the limit (5 — >■ after all calculations are performed. 
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As in Eq. (3.14), the general solution to Eq. (5.10) which satisfies the requisite boundary 
conditions at the stretched horizon reads 



V V{Vl)viV2)9rriVlPl \V2P2) =/ (Pout (Pl , ^^1 ) , Pout (P2 , ^^2) ) - /(Pout (P/l, ^^l) , Pout (P2, f 2) ) 

- /(Pout(Pl,^^l),Pout(P/l,f2)) + /(Pout(Pfc,^^l),Pout(Pft,^^2)) , 

(5.12) 

where as in the equilibrium case, f{xi,X2) is determined by the initial conditions specified 
in the distant past. 



To determine f{xi, X2) we invoke a similar argument used in Section III CI In particular, 
assuming initial data for grr{vipi\v2P2) was specified in the distant past, the only relevant 
initial data for evolution near v = comes from a very narrow strip which is very close to 
the event horizon. In this strip the relevant initial data is given by the coincident point limit 



of the symmetrized correlator as given in Eq. (3.15). Consequently f{xi,X2) must be given 
by ^ 

/(a;i,X2) = -— loglxi - X2I . (5.13) 
47r 



With 5'rr('^iPi|'^2P2) kuowu, Eq. ( |5.7 ) yields the symmetrized horizon correlator 



G%{vi,V2) 



TT 



(9^,,(9^2 log |1 - e 



(5.14) 



We will discuss the physical implications of this result in the next section. 

The horizon spectral density can be obtained along the same lines as the horizon 
symmetrized correlator with the subscript replacement rr ^ ra — ar in Eq. (5.7). Be- 



cause gra-ar{viPi\v2P2) satisfics the Same equations of motion and boundary conditions as 
grr{vipi\v2P2), the general solution for gra-ar{'Vipi\v2P2) takes the same form as Eq. (5.12). 



However, as in the equilibrium case discussed in Section |III C 2[ the initial conditions de- 
termining f{xi,X2) for the spectral density come purely from the equal time canonical 
commutation relations and yield 



/(xi, X2) = -sign(xi - X2) ■ 
With gra-ar kuowu, the horizon spectral correlator reads 



(5.15) 



Gra-ariVl,V2) = 2^/ri{vi)ri{v2)5' {vi - V2) 



(5.16) 



VI. SUMMARY 



Let us summarize our results. For definiteness we will describe how to compute string 
fiuctuations in the non-equilibrium geometry determined by Chesler and Yaffe We 
emphasize, however, that similar formulas can be used for other fields such as gravitons 
and dilatons. Fig. |6] shows the event horizon ro{v) together with the associated bifurcating 
outgoing geodesies of this geometry. The stretched horizon is located at r/j(f ) = ro{v) + e 

In contrast to a retarded Green function (which is a response to a source), a symmetrized 
correlation function, or a fiuctuation, is a time dependent correlation that arose from a 
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definite initial condition. In the case of Hawking radiation (see Fig. |2]), tliis initial condition 
is the result of an ultra-violet vacuum fluctuation (or a symmetrized correlation function) 
which originated close to the event horizon of the developing black hole. This UV fluctuation 
skims the event horizon following outgoing lightlike geodesies until late times. Then the 
fluctuation, which is no longer so ultra-violet, leaves the bifurcating horizon and induces 
stochastic motion in the string. We implemented this picture of Hawking radiation directly. 

We showed how the Hawking flux out of equilibrium can be packaged into an effective 
action on the stretched horizon which can be used to determine the effect of the Hawking 
radiation on the exterior dynamics. The horizon action through quadratic order is 



dvidv2 x'livi) [G':.,{vi,V2)\ x^'M , (6.1) 



where the horizon correlation functions are the horizon spectral density, 

G.Va.K,^2) = 2^^{vMv^)5'{v, - V,) , (6.2) 
and the horizon symmetrized correlator, 

Gt{vuV2) = - V^^(^^)^(^^) g,^g,^log|l-e-g^-"---^-')^-|. (6.3) 



vr 



Here x^{v) is the location of the string on the stretched horizon as a function of time. The 
coefficient r](v) = {\/X/2iT)gxxiv,rh{v)) determines the coupling between the world sheet 
fluctuations and the near horizon geometry, and the effective horizon temperature records 



the Lyapunov exponent of diverging geodesies along the bifurcating horizon (see Eq. (5.5) 



and Eq. (5.11)). More invariantly, it is related to the extrinsic curvature on the stretched 
horizon. In equilibrium, t]{v) is the drag of the heavy quark, and T^s is the Hawking 
temperature. The horizon spectral density (which is proportional to u) is determined by 
the canonical commutation relations of the 1 + lei effective theory which describes the near 
horizon dynamics. The horizon symmetrized correlator is determined by the initial density 
matrix of the effective theory far in the past. When the horizon action is coupled to the 
bulk, this action generates noise on the stretched horizon which induces the random motion 
of the quark in the dual field theory. For the stationary black hole the effective action for 
the string was determined previously using complicated analytic continuations pTj . 

The importance of these horizon functions is that they determine the spectral density 
and symmetrized correlations in the bulk and boundary theories. Indeed, the bulk spectral 
density (the commutator) and symmetrized correlator (the anti-commutator) are determined 
by propagating their horizon counterparts away from the stretched horizon 

iGra-ar{Viri\v2r2) = / dt;^ dt;2 G^a (^^l^^l l^^l^^/i (^^0 ) (^ra (t'2'^2 (^2) ) ^<^m-ar K 1^2) ) (6-4) 



Grriviri\v2r2) = J dv[dv2Graiviri\v[rhiv[)) Graiv2r2\v2rhiv2)) ^^^(^'11^2) • (6-5) 

Finally, these bulk correlation functions can be lifted to the boundary to determine the 
spectral density and symmetrized correlator in the field theory (see Appendix 0) . When the 
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fluctuations are thermalized, the two correlation functions satisfy the fluctuation dissipation 
theorem 



Thus we can numerically determine the fluctuations and monitor their approach to equi- 
librium using the formalism of this work. This numerical calculation will be presented in 
future work. 

Even without a complete numerical computation, some preliminary remarks can be made 
about equilibration in AdS^. The Wigner transforms of G^^(wi,f2) and G'fg^_^j.{vi,V2) obey 
the fluctuation dissipation theorem at high frequencyp] Specifically, for a typical non- 
equilibrium time scale r, we have 



where v = {vi + f2)/2. Thus the string is born into equilibrium at high frequency, and 
eventually frequencies of order the temperature and below equilibrate. This conclusion 
seems squarely aligned with the results of Ref. [29] which was limited to operators of high 
conformal dimension. However, it must be emphasized that the map between the stretched 
horizon and the boundary is non-trivial, especially for five-dimensional fields where the 
3-momentum influences the coupling between the field and the near horizon geometry. 

One popular picture of Hawking radiation is based on quantum tunneling (see, for in- 
stance, ^30j)- This picture relates the thermal factor in the emission rate to the change in the 
black hole entropy, which appears by way of the Euclidean action. It would be interesting 
to make contact with this picture using the effective action formalism. This would help 
make the universality of the result fully manifest. It would also be interesting to see if the 
simple non-equilibrium effective action presented in this paper finds a simple origin in the 
tunneling picture. 

The current derivation of Hawking radiation and correlation functions is similar to the 
2PI formalism of non-equilibrium field theory [31], and we hope this will make black hole 
physics accessible to a wider audience. The derivation uses the unstable nature of the bifur- 
cating event horizon to expand ultraviolet vacuum fluctuations. This exponential sensitivity 
to initial conditions has been called the transplanckian problem [32], and is characteristic 
of classically chaotic systems [33]. By exploiting this analogy, we hope that new insight can 
be found into the transplanckian problem and the Bekenstein entropy. Understanding the 
Bekenstein-Hawking entropy will require coupling the particle emission to the background 
metric, leading to a dynamical competition between quantum particles and the classical 
background field. The particle-field problem has been extensively studied in thermal field 
theory and the Color Glass Condensate [IHl El]- We hope to pursue these connections in 
future work. 

Comparison with recent literature: 



The Wigner transform is the Fourier transform with respect to the difference vi — V2, &s a function of the 
average, v = {vi + V2)/2. 




(6.6) 




(6.7) 
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Recently several papers appeared which addressed aspects of this paper as this work 
was being finalized [35]. First, a paper by Headrik and Ebrahim [35] used the equations 
of motion to solve for the symmetrized correlation function in AdS^-Va-idja. spacetimes. 
Headrik and Ebrahim reported on the "instantaneous thermalization" of AdSs. Second, a 
paper by Balasubramanian et al [22] computed the thermalization of operators with high 
conformal dimension by studying geodesies. Although the current paper is not limited to 
such operators, the basic conclusion that the field theory thermalizes first at high frequency 
is consistent with our conclusion about horizon Wigner transforms. However, as empha- 
sized above the map between the stretched horizon and the boundary involves physically 
important and non-trivial outgoing propagators. 
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Appendix A: Symmetrized correlations and the Keldysh formalism a lesson from 
the harmonic oscillator 

It is instructive in many respects to compute the symmetrized correlation function of the 
harmonic oscillator using the Keldysh formalism. The action of the oscillator is 

iSi — 182 = i dt [mXrXa — mulxrXa] , (Al) 

J to 



and is similar to the string action written in Eq. (4.4). In both systems we see that there 
are no rr type propagators in the action itself. This is because symmetrized correlations are 
the result of a correlation built into the initial state wave functions, i.e. the initial density 
matrix. 

In symmetrized type correlation functions the density matrix at an initial time to corre- 
lates the initial conditions for subsequent evolution. The density matrix for the harmonic 
oscillator in the ground state is '^o{xi)'^l{x2), and the symmetrized correlator is 

Grr{t, t)= I d<dx° ^o{x° + </2)^:« - x°/2) / Dx^Dx^e^^^'^^^x, (t)x, (t) . (A2) 

This is simplified by (i) introducing the Wigner transform, 

iy(x°,p°) = J dx°e-« ^oix: + x:/2)^:ix:-x:/2), 

(ii) integrating by parts in the action (with the boundary condition x^ — )• for t — )• +00), 

Jto 
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and finally (iii) integrating over all Xa, yielding 

r dx°dn° f 

Grr{t,t) = / W{x°,p°) 5{p" - mXr{to)) / DXr St [M Xr + mUoXr] Xr{t)Xr{t) . 

(A3) 

Here 5t[f{t)] denotes the functional delta function, S{f{t)). The meaning of this path 
integral is that Grr{t,t) is found by solving the equation of motion for a specified initial 
condition, 

po 

xAt) = x° cos(uJo(t - to)) H smiujoit - to)) , (A4) 

moOo 

and then averaging the square of this solution over the initial conditions specified by the 
Wigner transform. Performing this average for the ground state wave function of the oscil- 
lator reproduces the familiar result 

Grr{t, i) = I (0| {x{t),x{i)} |0) = cosMt - i)) . (A5) 

/ imujo 

The lesson from this analysis is that symmetrized correlation functions invariably arise from 
correlations in the initial density matrix which are propagated forward by the equations of 
motion. This dependence on the initial density matrix should be contrasted with retarded 
propagators which are independent of the wave function of the simple harmonic oscillator, 
i.e. 

9it - i) \xit),x{t)] = ~ ^ sm(uo(t - t)) , (A6) 

mtUo 

is a pure number. 



Appendix B: The Green Function Composition Rule 



In this appendix we detail the Green function composition rule stated in Eq. (3.6). In 
this section 1,2,3 denote the space-time points, e.g. 1 = (fi,ri). 

Suppose GraiM'^) and Gra(l|3) are retarded Green functions. Then the Wronskian of the 
two Green functions is 



W'^{2) = ^g^^Vhh>^'^{2) 



Gr,m—Gra{2\3) 



(Bl) 



^ y y ^ 

where d = d — d and is not intended to act outside of the square braces. Then, a short 
exercise shows that the divergence is 



d^W^{2) = aa(l|3)5(2, 3) - aa(3|l)5(2, V 



(B2) 



Assuming that 3 is inside the strip and 1 is outside the strip (see Fig. |3]) we can integrate 
over the strip to obtain the retarded Green function: 



aa(i|3) 



-§^W^{'2)= U'L,i^g^^Vhh^'^{2) 



GUl\2)—Gram 



(B3) 
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where dS^ is a surface surrounding the strip with outward directed normal, and the in- 
tegration is over space-time point 2. We next use the near horizon approximation for the 
leading factors (Eq. (2.6)), and neglect all surface terms except the integral over the stretched 
horizon. These surface integrals vanish because one of the Green functions vanishes. For 
instance, on the past surface (where v = —oo) (jra(2|3) must vanish since it represents the 
causal response at point 2 (past infinity) to a source at point 3. Since we have not specified 
the boundary conditions on the retarded Green function Gra, we are free to specify reflec- 
tive Dirichlet boundary conditions on the stretched horizon, i.e. take (jra(2|3) = Qrai'^l'^), 
as defined in the text. This specification does not interfere with the relevant initial data 
extremely close to the real event horizon. With this choice Eq. (B3) results in Eq. (3.6) 
given in the text. 



Appendix C: Prom bulk to bulk correlators, to boundary to bulk correlators, to 
CFT correlators 



Given a retarded bulk to bulk correlator GraiM'^') in the gravitational theory, we will 
show that the retarded correlator in the CFT is 



lim lim 

ri— >bnd r2— >bnd 



27T 



A 



2ti 



r2 



Graim 



(CI) 

In holography, the CFT correlation function is usually determined by setting the boundary 



conditions of five dimensional fields as r — )■ oo. We claim that Eq. (CI) is an equivalent 



prescription. Indeed, Eq. (CI) arises if the path integral discussion of the horizon effective 



action in Section |IVB| is applied to the boundary. Specifically, instead of integrating out 
the strip to find an effective action of the stretched horizon, one integrates out the entire 
bulk to find an effective action of the CFT. Put differently, from the viewpoint of the bulk 
theory, any small change in a boundary condition must be equivalent to inserting a suitable 
local operator close to that boundary, and for a Dirichlet condition that operator turns out 



to be a derivative of the field as in Eq. (CI) (see, for instance, [23]). 



To make direct contact with the pioneering work of Policastro, Son, and Starinets 
we define the boundary to bulk correlator f{tiri\t2) as a retarded solution which satisfies 



lim /(tirilta) = 5(ti-t2) 



J>bnd 



(C2) 



For static geometries in Fourier space this is usually called fuj{ri), and the retarded correlator 
in the CFT is usually 

+ ^g..^^h''dJUri)- (C3) 



CFT. 



CO 



2n 



Now we claim that the boundary to bulk propagator is simply related to the bulk to bulk 
propagator via 

fihn\t,)= hm ^Vhh'-{r2)dr,Graihn\t2r2). (C4) 



To show this we take ri fixed and large, and we integrate the equations of motions of the 



retarded Green function (Eq. (3.1)) with respect to the second argument over the pill-box 



25 



shown below. The radius of the lower surface is small compared to ri (but still large), while 
the radius of the upper surface is large compared to ri. 



bnd 



(C5) 



This yields 



lim lim ^-^—Qxx 

ri— >bnd r2— >bnd 27r 



v^/i^^(r2)T^G'(rA|r2t2) 
— lim lim ^^^g^x\fhK 

r2— >bnd ri— >bnd 27r 



(r2)^G(tiri|t2r2) =5(ti 
ar2 



(C6) 



The second term vanishes since GT.a^\'<'\\ii"<''i) obeys Dirichlet boundary conditions. Thus, 
the first term in Eq. (C6) obeys the same boundary conditions as the bulk to boundary 
propagator f{t\r'2\t2), Eq. (C2). Since both functions are retarded and obey the same 
equations of motion and boundary conditions, they are the same and Eq. (C4) holds. By 
extension, Eq. (CI) for the CFT correlators from bulk to bulk correlators is equivalent to 



the usual prescription. 
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